Effective temperatures of a driven, strongly anisotropic Brownian system 



O 

(N 

C 

00 
(N 



o 



i 

C 

o 
o 



> 

00 
00 

in 
in 

o 



X 



Min Zhang and Grzegorz SzamefB 
Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523, USA. 

(Dated: June 29, 2011) 

We use Brownian Dynamics computer simulations of a moderately dense colloidal system un- 
dergoing steady shear flow to investigate the uniqueness of the so-called effective temperature. We 
compare effective temperatures calculated from the fluctuation-dissipation ratios and from the linear 
response to a static, long wavelength, external perturbation along two directions: the shear gradient 
direction and the vorticity direction. At high shear rates, when the system is strongly anisotropic, 
the fluctuation-dissipation ratio derived effective temperatures are approximately wave- vector inde- 
pendent, but the temperatures along the gradient direction are somewhat higher than those along 
the vorticity direction. The temperatures derived from the static linear response show the same 
dependence on the direction as those derived from the fluctuation-dissipation ratio. However, the 
former and the latter temperatures are different. Our results suggest that the presently used for- 
mulae for effective temperatures may not be applicable for strongly anisotropic, driven systems. 

PACS numbers: 82.70.Dd, 05.70.Ln, 83.50.Ax 



I. INTRODUCTION 

Temperature is one of the most important parameters 
describing a macroscopic system in equilibrium. In con- 
trast, out of equilibrium, the very notion of temperature, 
its definition and its properties are still a subject of de- 
bate. This is true even for stationary nonequilibrium 
systems considered in this contribution. In principle, we 
can always define temperature in terms of average kinetic 
energy. It is not clear, however, whether such a definition 
leads to a parameter that has a more general significance. 
In addition, such a definition cannot be used for systems 
with overdamped stochastic dynamics {e.g. Brownian 
dynamics) (l[. 

In a seminal paper [2| Cugliandolo, Kurchan and Peliti 
proposed to define an effective, out-of-equilibrium tem- 
perature using a violation of a fluctuation-dissipation 
theorem. The starting point of their analysis was an ear- 
lier study [3| of the violation of the fluctuation-dissipation 
theorem in the low temperature phase of the spherical 
p-spin spin glass system. Building upon this study the 
authors of Ref. [2| showed that the ratio of the corre- 
lation and response functions, which in thermal equilib- 
rium is equal to the temperature, has some properties of 
temperature even in out-of equilibrium states. In partic- 
ular, they showed that an effective temperature defined 
through the ratio of the correlation and response func- 
tions controls the direction of the heat flow. 

Inspired by Ref. ||, Barrat, Berthier and Kurchan 
0-0 investigated the effective temperature defined in 
terms of the violation of the fluctuation-dissipation theo- 
rem for a classic glass-forming system, the so-called Kob- 
Andersen |7| binary Lennard- Jones mixture, undergoing 
a steady shear flow. A stationary drive has some tech- 
nical advantages compared to earlier studies of the effec- 



*Electronic address: szamel@lamar.colostate.edu 



tive temperature for the same system undergoing aging 
Q . The main one is that it can be expected that the sys- 
tem reaches a stationary non-equilibrium state in which 
time-translational invariance is restored. This fact makes 
easier both a computer simulation study and a theoreti- 
cal investigation. A very interesting feature of a system 
consisting of particles, as opposed to mean-field spin sys- 
tems, is that a great variety of response and correlation 
functions and the corresponding fluctuation-dissipation 
ratios can be defined. In principle, these different ra- 
tios could lead to a diverse set of effective temperatures. 
Then, it would be difficult to argue that any of these tem- 
peratures has a deeper meaning. However, Berthier and 
Barrat [6] showed that a class of fluctuation-dissipation 
ratios results in the same effective temperature. Further- 
more, they showed that the same temperature is regis- 
tered by a "thermometer", i.e. a subsystem coupled to 
the simulated system. 

In a related study Ono et al. @ investigated a model 
sheared athermal system (zero-temperature foam) . They 
considered a set of quantities that are equal to the tem- 
perature in equilibrium. Some of these quantities were 
based on static rather than dynamic linear response re- 
lations, and thus they were different from fluctuation- 
dissipation ratios considered in Refs. [2j-[6| from a funda- 
mental point of view. Ono et al. found that the quanti- 
ties they considered predicted the same value of effective 
temperature. Thus, these quantities could be used to 
define the effective temperature. 

In a follow-up study O'Hern et al. [10] addressed the 
question of the relation between effective temperatures 
defined on the basis of static linear response relations @ 
and the fluctuation-dissipation ratios 041] ■ They showed 
explicitly that in most cases these different definitions 
agreed over two and a half decades of effective temper- 
ature. However, in some cases they obtained different 
effective temperatures. The authors of Ref. [l(| noted 
that it is not entirely clear whether a given pair of conju- 
gate response and correlation function could be used to 
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define the effective temperature and which definition of 
the effective temperature (static or dynamic) should be 
used. 

We briefly mention here a com plet ely different system 
considered by Hayashi and Sasa 11]: a single Brownian 
particle moving in a tilted periodic potential. The advan- 
tage of this simple system is that many quantities can be 
calculated exactly [Tl|,[l2j]. In particular, the diffusion co- 
efficient and differential mobility can be calculated and 
their ratio can be used to define an effective tempera- 
ture in close analogy with fluctuation-dissipation ratios 
of Refs. [2|-l6|. Interestingly, Hayashi and Sasa showed 
that the effective temperature obtained from the ratio of 
the diffusion coefficient and differential mobility plays the 
role of the usual temperature in the large scale descrip- 
tion of their system. Consequently, they showed that, 
for this simple system, the fluctuation-dissipation ratio- 
derived temperature coincides with the temperature ob- 
tained from a long wavelength limit of a static linear 
response relation. 

Many-particle driven systems for which various defini- 
tions of effective temperatures investigated so far have 
two common features. The first one, emphasized for ex- 
ample in Refs. @, 4-6], is a separation of time scales: 
degrees of freedom evolving on short time scales seem to 
be thermalized with a different temperature than those 
evolving on the longest time scales. The latter degrees of 
freedom are thermalized with the effective temperature 
defined via the fluctuation-dissipation ratio. It should be 
noted that the importance of the time scale separation 
is somewhat less clear for the static definitions of the ef- 
fective temperature. The second common feature is that 
driven systems considered in Refs. 0-SSE3] were struc- 
turally isotropic. In particular, shear rates used in Ref. 
@ were small enough for the steady state static structure 
factors in the flow, velocity gradient and vorticity di- 
rections to be identical within simulational uncertainty. 
Thus, there was no reason to expect that the effective 
temperature defined using wave-vector-dependent corre- 
lation functions would depend on the direction of the 
wave- vector and thus be anisotropic. 

Here we consider a moderately dense model colloidal 
system undergoing a shear flow in the limit of high shear 
rates. As a result of the high shear rate our system is 
strongly anisotropic. In particular, the static structure 
factor in the velocity gradient direction is considerably 
larger than its equilibrium value [l|| although no tran- 
sition to a layered state [l5[ is observed. We test three 
definitions of effective temperature along two different 
directions, the velocity gradient direction and the vor- 
ticity direction. The first two definitions are based on 
the fluctuation-dissipation ratios. Specifically, the first 
definition is based on a violation of the Einstein rela- 
tion between the friction coefficient and the self-diffusion 
coefficient, and the second definition is based on the vio- 
lation of the fluctuation-dissipation relation between the 
transient response describing the change of the tagged 
particle density due to an external potential acting on 



the tagged particle and the self-intermediate scattering 
function. The third definition is a static one based on 
the violation of the linear response relation describing 
the change of the density profile due to a static, long 
wavelength, external potential. 

All three definitions allow us to calculate effective tem- 
peratures pertaining to different directions relative to the 
shear flow. In addition, the second definition allows us to 
calculate effective temperatures corresponding to differ- 
ent wave- vectors. We find that the effective temperatures 
derived from the fluctuation-dissipation ratios are wave- 
vector independent, but at high shear rates the tempera- 
tures along the velocity gradient direction are somewhat 
larger than those along the vorticity direction. 

The dependence of the effective temperatures on the 
direction correlates with the fact that at high shear 
rates the system is strongly anisotropic. In addition, we 
find that the effective temperatures obtained from the 
static linear response relation show the same dependence 
on the direction as those derived from the fluctuation- 
dissipation ratios. However, the effective temperatures 
obtained from the static linear response relation are dif- 
ferent from those derived from the fluctuation-dissipation 
ratios. Finally, we show that if the static linear response 
definition is generalized to shorter wavelengths (i.e. to 
more rapidly varying in space periodic potentials), the 
resulting effective temperatures have complicated wave- 
vector and direction dependence. 

The paper is organized as follows. In the next sec- 
tion, Sec. [H] we describe the details of the simulation. 
We discuss the structure of the sheared system in Sec. 
IIII1 in particular, we show that at high shear rates the 
system is strongly anisotropic. In Sees. IIV1 IVl and [Vl] 
we present effective temperatures determined using the 
three different definitions. We close the paper with a 
short discussion presented in Sec. I VIII In the appendix 
we examine the wave- vector dependence of the effective 
temperatures defined through the static linear response 
relation. 



II. 



SIMULATION DETAILS 



We simulated one of the systems investigated in an 
earlier publication ]16]. Briefly, the system consists of 
N = 1372 spherical colloidal particles. We used the 
screened Coulomb potential originally introduced by Ras- 
togi et al [l5j : 



V{r) = A 



exp (—k(t — a)) 



(1) 



with A = A7hksTa and k<j — 24. The cutoff distance 
was 3cr. We simulated the system at dimensionless den- 
sities, (N/V)cr 3 = ner 3 equal to 0.408 (V is the system's 
volume) . We previously established [l6| that this density 
corresponds to an effective hard-sphere volume fraction 
equal to 0.43. Thus, our system is moderately dense. 



We performed Brownian dynamics simulations. The 
equation of of motion for the ith particle is, 
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l 5ri £y(| ri _ r .|) +7?i ( t ) (2) 



where 7 is the shear rate, yi is the coordinate of parti- 
cle i along y direction, e x is the unit vector along the x 
direction, £ is the friction coefficient of an isolated par- 
ticle, and the random noise rji satisfies the fluctuation- 
dissipation theorem, 



(Vi(t)r] j {t'))=2DoS(t-t% j I. 



(3) 



In equation ([3]), I is the unit tensor and the diffusion 
coefficient Dq = fc^T/^o where ks is Boltzmann's con- 
stant. We will present the results in terms of reduced 
units with a and a 1 / Dq being the units of length and 
time, respectively. 

The equations of motion, Eq. ([2]), were solved with 
Lees-Edwards boundary conditions for shear rates 7 = 2, 
5, 10 and 20. To solve the equations of motion we used 
a second order Brownian dynamics algorithm [17| | with a 
time step of 10~ 4 for shear rates 2 and 5, and a time step 
of 5 x 10 -5 for shear rates 10 and 20. 



III. 



STRUCTURE OF THE SHEARED SYSTEM 



A fluid's structure is usually discussed in terms of its 
pair distribution function. Under shear, the pair distri- 
bution function <?(r), 



9(r) 




(4) 



is anisotropic (hereafter the brackets (...) denote the non- 
equilibrium steady state average). The anisotropic pair 
distribution is somewhat difficult to visualize and usu- 
ally only its two lowest order projections on spherical 
harmonics are monitored llal: 



3 s 0) = J <%(r), g + (r) = ^ J drxyg(v). 



(5) 



Here we follow notation introduced by Hess and Hanley 
19]: g s (r) denotes the isotropic component of the pair 
distribution function, which in the absence of the shear 
flow reduces itself to the familiar equilibrium pair distri- 
bution function g eq (r)- Next, g+(r) denotes the projec- 
tion of the pair distribution function onto xy. In princi- 
ple, the anisotropic pair distribution g(r) can be recon- 
structed from the projections 



g(r) = g s {r)+g + (r)xy + 



(6) 



but, in practice, for high shear rates the expansion in- 
dicated above is slowly convergent and the higher order 
terms, indicated by dots in Eq. ([5]), are important (20j . 
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FIG. 1: (Color online) The isotropic component of the pair 
distribution function, g 3 (r), as a function of the particle 
separation, r, for different shear rates. 



In the limit of low shear rates g s = g eq 4- 0(7), g+ is 
proportional to 7 and all projections other than g s and 
g + are 0(7). In addition to being the dominant small 7 
correction to the equilibrium pair distribution, the im- 
portance of g+ comes from the fact that for pair-wise 
additive interactions it determines the xy component of 
the interaction contribution to the stress tensor, <j xy , 



'xy 



n 2 47T 
~2~15 



dr 



9+(r) 



(7) 



and, as a consequence, the non-linear viscosity 77(7) = 

We show g s and g+ in Figs. Q] and [2] With increas- 
ing shear rate the isotropic component, g s , changes lit- 
tle whereas the amplitude of g+ systematically increases. 
The inset in Fig. [5] shows that this increase of g+ is sub- 
linear in 7. This fact is a reflection of the well-known 
shear thinning phenomenon: sub-linear increase of gf+ 
leads to a sub-linear increase of the xy component of the 
stress tensor and a decrease of the non-linear viscosity 
with increasing shear rate. 

While the isotropic component of the pair correla- 
tion function g s changes little with increasing shear, at 
the largest shear rate investigated, 7 = 20, the first 
anisotropic component g + is comparable to g s . In fact, 
although it is not evident from Fig. [5J the system is 
strongly anisotropic already at 7 = 5. The first evidence 
of this anisotropy can be obtained by investigating the 
static structure factor S(k), 

1 N 
5(k) = - (p(k)p(-k)> , p(k) = £ expHk • r 4 ), (8) 

i=l 

along specific directions in k space. 

In Fig. [3] we show the structure factor along two dif- 
ferent directions: the velocity gradient direction e y and 
the vorticity direction e z . We see that the first peak 
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FIG. 2: (Color online) The first anisotropic component of the 
pair distribution function, g+(r), as a function of the parti- 
cle separation, r, for different shear rates. The inset shows 

g+(r)/i- 
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FIG. 3: (Color online) The structure factor, S(k), for k along 
two different directions for different shear rates, (a) gradient 
direction, e y . (b) vorticity direction, e z . 



of the structure factor along the gradient direction in- 
creases strongly with shear rate starting at a shear be- 
tween 7 = 2 and 7 = 5. In contrast, the structure factor 
along the vorticity direction shows almost no shear rate 
dependence. Similar behavior was seen before for a sys- 
tem under oscillatory shear flow [l4[ . The increase of the 
structure factor along the velocity gradient direction can 
be interpreted in terms of local ordering into a layered- 
like structure fll| . 

To investigate the cross-over between strong shear rate 
dependence of the structure factor along the gradient di- 
rection and weak shear rate dependence along the vortic- 
ity direction it is useful to visualize the structure factor 
in some planes in the wave-vector space. We investi- 
gated the shear rate dependence of the structure factor 
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FIG. 4: (Color online) The structure factor, S'(k), in the 
velocity-gradient plane, for different shear rates. 
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FIG. 5: (Color online) The structure factor, S(k), in the 
gradient-vorticity plane, for different shear rates. 



in two planes, the flow-velocity gradient plane (x-y) and 
the gradient-vorticity (y-z) plane. Fig. [4] shows that 
the increase of the structure factor along the velocity 
gradient direction is quite localized in the flow-velocity 
gradient plane. Moreover, we note that at the highest 
shear rate the whole ring-like pattern becomes strongly 
deformed. In contrast, Fig. [3] shows that in the veloc- 
ity gradient-vorticity plane the transition from strong in- 
crease of the structure along the gradient direction to 
almost shear rate independent structure factor along the 
vorticity direction is more gradual. In addition, in the 
velocity gradient-vorticity plane the ring-like structure is 
still quite un-deformed even at the highest shear rate. 

To summarize, in this section we showed that at 
high shear rates the steady-state system is strongly 
anisotropic. 



IV. EFFECTIVE TEMPERATURE: EINSTEIN 
RELATION 

According to the celebrated Einstein relation, the 
mean-square displacement of the tagged particle is re- 
lated to the tagged particle motion under the influence 
of a weak external potential acting on the tagged particle 
only. Specifically, for a system in equilibrium, in the long 



time limit we can define the equilibrium self-diffusion co- 
efficient D via the tagged particle mean square displace- 
ment, 

lim l/| ri (t)- ri (0)| 2 ) = 6D eq , (9) 

t— >oo t \ I eq 

where the brackets (•••} e „ denote the equilibrium aver- 
age. In addition, we can define the equilibrium self- 
mobility coefficient fi via the systematic displacement of 
the tagged particle under the influence of the weak ex- 
ternal force of magnitude Fq acting along an arbitrary 
direction specified by the unit vector f , 

Pb-s-0 t-too tt<o \ I eq 

where {■■■)^ q denotes an average in a system without 
shear, subjected to a constant force acting on the tagged 
particle. 

D eq and /i eq are connected via the so-called Einstein 
relation involving the system's temperature T, 

D eq = Tfleg. (II) 

For a steady state system under shear both the self- 
diffusion coefficient and the self-mobility coefficient de- 
pend on the shear rate. We note here that, for col- 
loidal systems under oscillatory shear, the shear rate de- 
pendence of the former was studied over two decades 
ago in experiments [21] and computer simulations [Tij . 
These studies were followed by theoretical studies of self- 
diffusion in semi-dilute [22j and concentrated [23| col- 
loidal systems under steady shear. In the former work, 
the self-mobility was also studied and the violation of the 
Einstein relation (TTT|) was noted. The realization that 
Eq. (fTTj) is violated was not followed, however, by the 
crucial step of the identification of an effective temper- 
ature. Berthier and Barrat [H, H3] made this last step, 
turned Eq. (fTTj) around, defined and calculated the ef- 
fective temperature 

Teff = D/p. (12) 

It should be emphasized that both quantities at the right- 
hand-side of Eq. (|12l) depend on the shear rate j but in 
a slightly different way. It is this last fact that gives rise 
to a non-trivial effective temperature. 

Since our system is anisotropic, we investigated two 
different effective temperatures based on the violation 
of Eq. (fTT]) along two different directions, the gradient 
direction and the vorticity direction, 

~^eJ7 = Dyy/fJ-yy, (13) 
Teff = D zz /flzz- (14) 

Here D yy and D zz are the diagonal components of the 
self-diffusion tensor along the gradient and vorticity di- 
rections, respectively, 

lim - ((oi(t) - a t {0)f) = 2D aa , (15) 




5 10 15 20 5 10 15 20 



Y 

FIG. 6: (Color online) (a-b): shear rate dependence of the 
diagonal components of the self-diffusion and mobility ten- 
sors. Note that in equilibrium, i.e for 7 = 0, for our value of 
temperature, T = 1, the self-diffusion and the mobility coef- 
ficients are equal, D eq — fi eq . (c-d): shear rate dependence of 
the effective temperature obtained from the self-diffusion and 
mobility tensors, (a) and (c): gradient direction, e,,. (b) and 
(d): vorticity direction, e z . 

where a = y or a = z. Furthermore, in Eqs. (I13H14[) fiyy 
and n zz are the diagonal elements of the self-mobility 
tensor 

lim lim — ((a t (t) - a l (0))) F ° = ^ aa , (16) 

F ->0 i->oo tl<o 

where, again, a = y or a — z, and the weak external 
force of magnitude Fq is applied on the tagged particle 
in the direction a. 

The elements of the self-diffusion tensor can be easily 
calculated. In contrast, an evaluation of the elements of 
the self-mobility tensor is a little more challenging. We 
followed the procedure described in Ref. [6J and intro- 
duced a constant force acting on 686 randomly selected 
particles of the system, with the force in the positive 
direction along the relevant axis for 343 particles and 
the force in the negative direction for the remaining 343 
particles. The average systematic displacement of these 
particles, weighted by the sign of the force, is a good 
approximation for ((a,(f) — o,(0))) and can be used to 
calculate fi aa . We checked this assertion by re-doing this 
calculation with a constant force acting on 340 randomly 
selected particles of the system. As expected, the statis- 
tics was worse, but the results were compatible with the 
calculation using a constant force acting on 686 particles. 

The results are presented in Fig. [BJ The upper panels 
of this figure show that both the diagonal elements of the 
self-diffusion tensor and the diagonal elements of the self- 
mobility tensor increase with increasing shear rate. Both 
quantities are somewhat larger along the gradient direc- 
tion than along the vorticity direction. We note that for 
a system under an oscillatory shear an early simulation of 
Xue and Grest [13] found the opposite behavior for the 
elements of the self-diffusion tensor. In contrast, for a 
system under steady shear a theoretical study of Indrani 
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and Ramaswamy [23j found that the diagonal elements 
of the self-diffusion tensor along the gradient direction 
was larger than those along the vorticity direction. This 
qualitatively agrees with the results for the elements of 
the self-diffusion tensor showed in Fig. [6j It should be 
emphasized, however, that in Ref. [23] the violation of 
the Einstein relation was neglected and thus Indrani and 
Ramaswamy's approach cannot describe the difference 
between self-diffusion and self-mobility tensors. 

The lower panels of Fig. [6] show that, as found in 
earlier investigations, the effective temperature increases 
with increasing shear rate. The temperatures along the 
gradient and vorticity directions are quite close, but the 
one along the gradient direction is systematically some- 
what larger than the one along the vorticity direction. 
We should note, however, that the uncertainty of these 
effective temperatures prevents us from making a defini- 
tive statement that the two temperatures are different. 



V. EFFECTIVE TEMPERATURE: FDR 
INVOLVING THE SELF-INTERMEDIATE 
SCATTERING FUNCTION 

According to the standard fluctuation-dissipation rela- 
tion, in equilibrium the susceptibility describing a time- 
delayed change of the tagged particle density is propor- 
tional to the self-intermediate scattering function, 



Xs, eq (k:t) = i[l-F. ia ,(fc;t)] 



(17) 



Here Xs,eq(k; t) is the susceptibility describing the change 
of the Fourier component k of the tagged particle density 
due to a weak external force acting on the tagged particle, 



d 

Xs.eqik; t)=gr (exp [ik ■ Ti(t)]) ei 



q \h=o 



(18) 



where (■■■) eq denotes an average calculated in a system 
that was in equilibrium at t = and subsequently was 
subjected to a constant in time force Ff xt acting on the 
tagged particle only, 

Ff* = -d rt [2fc cos (k •!■<)] . 

Furthermore, F se q{k;t) in Eq. (|17[) i s the equilibrium 
self-intermediate scattering function, 

F s , eq (k;t) = (exp[ik- (r,(t) -^(0))]). (19) 

For a system undergoing a steady shear flow, the rela- 
tion (fl~T|) is violated. Following @ we define an effective 
temperature through the violation of this relation, 



I eff 



1-F s (k a ;t) 
Xs(k a ;t) 



(20) 



All quantities at the right-hand-side of Eq. (|2"0)l are cal- 
culated for a system undergoing a steady shear. Specifi- 



cally, F s (k a ;t) is the steady state self-intermediate scat- 
tering function. Moreover, Xs(k a ]t) describes the time- 
delayed change of the tagged particle density, 



Xs{k;t) 



d_ 

dh 



[exp[ik ■ ri(t)]) 



(21) 



where (...) denotes an average calculated in a system 
that was undergoing a steady shear at t — 0, and subse- 
quently was subjected to a constant in time force acting 
on the tagged particle. It should be noted that the ef- 
fective temperatures in definition (|20[) , in addition to the 
dependence on the direction of the force, a = y, z, can 
also depend on the magnitude of the wave- vector (to sim- 
plify notation this possible dependence is not indicated 
explicitly). As discussed in the Introduction, the authors 
of Ref. [6] found that, for a supercooled system under- 
weak shear, the effective temperatures defined through 
Eq. (I2TJ1) does not depend on the wave- vector. 

The calculation of the self-intermediate scattering 
function is straightforward. In contrast, it is more dif- 
ficult to calculate the susceptibility. We followed the 
procedure described in Ref. @. Specifically, we con- 
sidered a system with a constant force imposed on 686 
randomly selected particles of the system, with the force 
in the positive direction along the relevant axis for 343 
particles and the force in the negative direction for the 
remaining 343 particles. Furthermore, in order to evalu- 
ate the susceptibility using the un-perturbed trajectories, 
we followed a recently introduced approach [25j and we 
performed the derivative in Eq. (|21[) before taking the 
average. This necessitates solving additional equations 
of motion for quantities dvijdh [25|], but it allowed us to 
calculate the susceptibility from one long, unperturbed 
trajectory. It should be noted that, in contrast to the 
calculation described in Ref. [25| . we did not have to 
perform several independent runs. This simplification 
follows from the fact that while in Ref. [25| a response 
of an aging system was studied, we are interested in a 
response of a steady state system. 

In Fig. [7]we show parametric plots of the susceptibility 
versus the correlation function for several wave- vectors 
along the y and z axes at the highest shear rate, 7 = 20. 
In plots of this type the effective temperature is the (mi- 
nus) of the inverse slope of the data. Moreover, short 
times correspond to the lower right corner and long times 
to the upper left corner of the plot. One can see from 
Fig. [7]that at short times the fluctuation-dissipation (Ti~7|) 
is obeyed. In contrast, at long times an effective temper- 
ature higher than the heat bath temperature is observed. 
Notice that at long times, for a given direction, the slopes 
of the lines fitted to the data obtained at different mag- 
nitudes of the wave-vectors seem very close. This agrees 
with the findings of Ref. [6j. In contrast, the slopes of 
the lines fitted to the data obtained at two different di- 
rections are different. 

The results of the analysis of parametric plots at all 
shear rates investigated are presented in Fig. [5] In 
agreement with the behavior found from the analysis of 





C(t) 

FIG. 7: (Color online) Parametric plots \ s (k a ;t) vs. F 3 (k a ;t) 
for wave-vectors along two directions, (a) gradient direction, 
e y . (b) vorticity direction, e z . The magnitudes of the wave- 
vectors are k a = 4.61 5.45, 7.55 and 10.48, listed from bottom 
to top. The solid lines have slopes equal to — 1/T = — 1. The 
dashed lines have slopes equal to — 1 /T" s j . The insets indicate 
the magnitudes of the wave-vectors in relation to the steady 
state structure factors for wave- vectors along the gradient and 
vorticity directions. Recall that at this high shear rate the 
steady state structure factors along these directions are quite 
different, see Fig. [3] 



the violation of the Einstein relation, we find that that 
in both directions and at all wave-vectors the effective 
temperatures increase with increasing shear rate. The 
temperatures obtained from different magnitudes of the 
wave-vectors are reasonably consistent. This agrees with 
the finding of Berthier and Barrat @. As found in the 
previous section, temperatures obtained from the exter- 
nal force along the gradient and vorticity directions are 
close, but the one along the gradient direction is system- 
atically larger than the one along the vorticity direction. 



FIG. 8: (Color online) Shear rate dependence of effective 
temperatures obtained from the violation of the fluctuation- 
dissipation relation involving the self-intermediate scattering 
function, (a): gradient direction, e y . (b): vorticity direction, 
e z . 



VI. EFFECTIVE TEMPERATURE: STATIC 
LINEAR RESPONSE DEFINITION 

For a system in equilibrium, the static structure factor 
plays the role of the response function: the change of the 
average density due to a weak, static, periodic in space, 
external potential, 

V(r) = V sm(q-r), (22) 

which has been turned on at t = — oo, is given by the 
following well known relation: 

5n(r) = -^My sin(q.r) (23) 

where T is the system's temperature and S eq (q) is the 
equilibrium structure factor. 

The main finding of Ref. [l(| is that the effective 
temperature obtained from the long wavelength (i.e. 
small wave-vector) limit of some linear response rela- 
tions agrees with that obtained from the fluctuation- 
dissipation ratios. Inspired by this result here we use the 
relation (f2"3"|) to define an effective temperature. Specif- 
ically, we assume a sheared, stationary suspension, im- 
pose a slowly varying external potential (f22j) (we use the 
smallest wave- vector allowed by periodic boundary condi- 
tions), wait a time long enough for the system to reach a 
new stationary state, and measure the density change due 
to the external potential. We use two different strengths 
of the effective potential (i.e. two different values of Vo) 
to extrapolate to the linear response limit implicit in Eq. 
(|23|) . In the limit of weak external potential the density 
change is proportional to Vo, 

5n(r) = — na(q)Vo sin(q • r) (24) 

where we made explicit the fact that under shear the 
density change depends on the orientation of the vector q 
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FIG. 9: (Color online) (a-b): shear rate dependence of the 
steady state structure factor S(q a ) and the response coef- 
ficient a(q a ), a — y,z, where q a is the smallest non-zero 
wave-vector allowed by the periodic boundary conditions, 
q y = qz — 0.42. (c-d): shear rate dependence of the effective 
temperature obtained from the static linear response, (a) and 
(c): gradient direction, e y . (b) and (d): vorticity direction, 



with respect to the shear flow. Next, we use the following 
generalization of the relation ([23]) to extract the static 
effective temperature: 



T a = S(q a ) 
Cff a(q a ) 



(25) 



It should be noted that the definition (|2"5)) involves the 
shear rate dependent, anisotropic structure factor S(q). 
This agrees with the static definition adopted in Ref . ■ 
Thus, the static effective temperature (|25|) is a ratio of 
two shear rate-dependent quantities. 

In contrast to the procedure adopted in Ref. [l(| we 
did not monitor the time dependent (transient) density 
change after the external potential was turned on, and 
thus do not obtain the static effective temperature from 
an intercept of a parametric plot of the response (i.e. 
the density change) as a function of the time dependent 
density correlation function. However, our definition of 
the effective temperature is equivalent to that used in 
Ref. ftp. 

In Fig. [5] we present the shear rate dependence of 
the structure factor and the linear response coefficient a. 
Since we are interested in the response to a long wave- 
length perturbation, both quantities are shown for the 
smallest wave-vector compatible with periodic boundary 
conditions, q y = q z = 0.42. While both quantities de- 
crease with increasing shear rate, the coefficient a has 
a somewhat stronger dependence on j. As a result, as 
shown in lower panels in Fig. [51 the effective tempera- 
tures calculated for the perturbation along the gradient 
and vorticity directions increase with shear rate. The ef- 
fective temperature along the gradient direction is larger 
than that along the vorticity direction. This agrees with 
the behavior obtained from the Einstein relation and the 




Ml: Einstein Relation (M 
M2: FDR ' ' 

M3: linear response 
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FIG. 10: (Color online) Comparison of the shear rate de- 
pendence of effective temperatures obtained using three dif- 
ferent definitions. Dotted line: Einstein relation; solid line: 
violation of the fluctuation-dissipation relation for the self- 
intermediate scattering function at the wave-vector corre- 
sponding to the peak position of the structure factor; dashed 
line: static linear response to a long wavelength perturbation, 
(a) gradient direction, e y . (b) vorticity direction, e z . 



fluctuation-dissipation ratios. However, the values of the 
effective temperature obtained from static linear response 
are different than those obtained from the fluctuation- 
dissipation ratios. 



VII. DISCUSSION 

To investigate the uniqueness of the notion of the so- 
called effective temperature of a driven colloidal system 
we studied three definitions of such a temperature. The 
main results are summarized in Fig. 1101 

We evaluated the effective temperature obtained from 
the violation of the Einstein relation Q and the ef- 
fective temperature obtained from the violation of 
the fluctuation-dissipation relation involving the self- 
intermediate scattering function [6| . The latter definition 
allows to calculate the effective temperature at different 
wave- vectors. Formally, the former definition can be con- 
sidered as the limiting case of the latter one for vanishing 
magnitude of the wave-vector. Berthier et al. [6| found 
that these two definitions resulted in compatible values 
of the effective temperatures. Our results suggest that 
at higher shear rates, when the system is anisotropic, the 
two definitions give the same result along the gradient 
direction and along the vorticity direction. However, the 
effective temperatures along the gradient direction are 
somewhat larger than those along the vorticity direction. 

In addition, we investigated a definition of the effec- 
tive temperature based on a static linear response rela- 
tion. The motivation for this was two-fold. First, sim- 
ilar relations were investigated by O'Hern et al. fToj j . 
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Second, Hayashi and Sasa found that for their sim- 
ple, single particle system, the effective temperature de- 
fined through a fluctuation-dissipation ratio determines 
the response of the system to a long wavelength potential. 
We found that the effective temperatures obtained from 
the static linear response relation were different from 
those obtained from violation of the Einstein relation and 
of the fluctuation-dissipation relation involving the self- 
intermediate scattering function. However, in the long 
wavelength limit static linear response derived effective 
temperatures had the same dependence on the direction: 
temperatures along the gradient direction were somewhat 
larger than those along the vorticity direction. In the ap- 
pendix we extended the static linear response definition 
of the effective temperature to shorter wavelength per- 
turbations. As shown in Fig. [T^l at shorter wavelengths 
(larger wave-vectors) the static linear response relation 
leads to a variety of different effective temperatures, and 
in some cases with dependence on the direction opposite 
to that shown in Fig. [TU1 

It would be interesting to repeat our numerical inves- 
tigation for a more strongly interacting system with a 
clear separation of time scales. To this end, however, one 
would have to use a mixture (possibly, Kob- Andersen bi- 
nary Lennard- Jones mixture) because we found that our 
single component system forms layers at high density and 
large shear rates (as was first reported in Ref. (HI). 

It would also be interesting to develop some theoretical 
understanding of the definitions of the effective tempera- 
ture investigated here. First steps in this direction have 
already been taken [24| [26|. It is still not clear, how- 
ever, when one has to use a fluctuation-dissipation based 
definition of the effective temperature and when one can 
use a static linear response based definition for a many- 
particle, strongly interacting, sheared system. 

More generally, it would be worthwhile to establish 
what, if any, properties of a many-particle, strongly in- 
teracting, sheared system are influenced by which defi- 
nition of the effective temperature. Again, first steps in 
this direction have already been taken [23|. We would 
like to mention here a very interesting result concerning 
driven systems with Newtonian dynamics: alignment of 
particles in a shear flow depends on whether kinetic or 
configurational [lj temperature is thermostated (28| . 

The main conclusion from our study is that for a 
strongly driven, anisotropic system presently used, defini- 
tions of the effective temperature may not give a unique, 
direction-independent result. It is, therefore, imperative 
that the realms of validity and importance are established 
for these definitions. Only then we will be able to fully 
assess the relevance of the existing definitions. 
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FIG. 11: (Color online) (a-b) and (e-f): shear rate depen- 
dence of the structure factor, (c-d) and (g-h): shear rate 
dependence of the linear response coefficient, (a), (c), (e) and 
(g): gradient direction, e^. (b), (d), (f) and (h): vorticity 
direction, e z . (a-d) show smaller wave- vectors and (e-h) show 
larger wave-vectors. 



Appendix 

Here we generalize the static linear response based def- 
inition of effective temperature. The main motivation for 
this generalization is that, as shown in Fig. [TU1 the static 
linear response to a long wavelength perturbation gives 
effective temperatures that are different from those ob- 
tained using the fluctuation-dissipation ratios. This dif- 
ference prompted us to check whether it is still present 
at larger wave- vectors. 

Effective temperatures presented in this section were 
obtained using relation (f25j) generalized to larger wave- 
vectors: 
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In Eq. (|26|) the effective temperature may depend on the 
wave- vector. 
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FIG. 12: (Color online) Shear rate dependence of the static 
effective temperature at a number of wave-vectors, (a) and 
(c) : gradient direction, e y . (b) and (d): vorticity direction, 
e z . (a-b) show smaller wave- vectors and (c-d) show larger 
wave- vectors. 



In Figs. [TT] and [12] we show the shear rate dependence 
of the structure factor, linear response coefficient, and ef- 
fective temperature for different wave- vectors along two 
directions, the gradient direction and the vorticity di- 
rection. We find that effective temperature as defined 
through Eq. (1261) depends significantly on the wave- 
vector. In particular, whereas at small wave- vectors 
temperatures along the gradient direction are somewhat 
larger than those along the vorticity direction, at some 
of the larger wave-vectors this ordering is reversed. Our 
finding emphasizes the need for some theoretical under- 
standing of the static linear response definition of effec- 
tive temperature. 
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